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The proton and neutron inclusive distributions in the projectile fragmentation region of pA colli- 
sions are studied in the valon model. Momentum degradation and flavor changes due to the nuclear 
{3J[)' medium are described at the valon level using two parameters. Particle production is treated by 

^ , means of the recombination subprocess. The centrality dependences of the net proton and neu- 

■ tron spectra of the NA49 data are satisfactorily reproduced. The effective degradation length is 

' determined to be 17 fm. Pion inclusive distributions can be calculated without any adjustable 

parameters. 
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I. INTRODUCTION 



> 

' The study of proton-nucleus (pA) collisions is important because they are tractable intermediaries between pp and 
, AA collisions, when intense interest exists in discovering the extent to which the dense medium created in an AA 
collision differ from that of linear superpositions of pp collisions. One of the properties of pA collisions is that the 
momentum of the leading baryon in the projectile fragmentation region is degraded, a phenomenon commonly referred 
to, somewhat inappropriately, as baryon stopping. Such a transference of baryon number from the fragmentation to 
1-^, I the central region contributes to the increase of matter density at mid-rapidity, thereby raising the likelihood of the 
I . formation of quark-gluon plasma. Thus it is important to understand the process of baryon momentum degradation 
^ ■ and its dependence on nuclear size or centrality. 

^ , Since it is not feasible to perform first-principle parameter-free calculations of the momentum degradation at this 
' point, experimental guidance is of crucial importance. Recently, several experiments have produced useful data on 
the subject, in particular, E910 and E941 at the AGS, and NA49 at the SPS 0. It is the xp dependences of the 
distributions of p — p and n — n that we shall focus on; moreover, their dependences on centrality will guide us in our 
determination of the nuclear effect on baryon momentum degradation. 

The degradation of baryon momentum has been studied in various approaches before In Refs. the 

investigations are done at the nucleon level, while in Refs. the string model is the basis. Since it is questionable 

that the concept of color strings can be relevant in heavy-ion collisions where the abundance of color charges in the 
overlap region renders unlikely the development of constricted color flux tubes , our approach in this paper will be 
on the various levels of the constituents of the nucleon that are consistent with the parton model. More specifically, 
we shall use the valon model ||T^-||l^ to keep track of the momenta of the constituents and the recombination model 
[T^ to describe the hadronization of the partons. The valons play a role in the collision problem as the constituent 
quarks do in the bound-state problem. Thus a nucleon has three valons which carry all the momentum of the 
nucleon, while each valon has one valence quark and its own sea quarks and gluons. Although soft processes are 
non-perturbative, the valon model (including recombination) nevertheless provides a systematic way of calculating all 
subprocesses that contribute to a particular inclusive process. Some of the subprocesses can be identified with certain 
diagrams in other approaches, e.g., baryon junction and diquark breaking terms |@, ^. 

The valon distributions in the proton will be determined by fitting the parton distributions at low . The effect of 
the nuclear medium on the valon distribution of the proton projectile will involve two parameters, one characterizing 
the momentum degradation and the other flavor flipping. The color indices are all averaged over, since multiple 
gluon interactions, as the projectile traverses the target nucleus, are numerous and uncomputable. Their effects are, 
however, quantified in terms of the number, of target nucleons that participate in the pA collision. For that reason 
it is important that the experimental data must have centrality selection expressible in terms of the average v. 

The inclusive process p + A ^ p + X has been studied before in the valon model . However, the pertinent data 
available at that time had no centrality cuts and were for fixed pT fl^ . Now, NA49 has p-Pb data on the production 
of p, p, n and n in the proton fragmentation region (with contribution from the nuclear target subtracted) for various 
values of P Q . To extract more information from these refined data we have improved our formulation of the valon 
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model, allowing for momentum degradation of each of the valons independently and for changes in their flavors. The 
effect on the nucleon momentum distributions can be calculated and the effective degradation length deduced. 

In Sec. II we determine the valon and parton distribution functions from the published distribution functions 
at low . The valon model for pA collisions is then discussed in Sec. Ill, in which momentum degradation is 
formulated. In Sec. IV we consider in detail the projectile fragmentation process and the subsequent process of 
quark recombination so that the inclusive distribution of nucleons can be calculated. The net proton and neutron 
distributions are then calculated and compared to the experimental data in Sec. V. Predictions for pion production 
in the proton fragmentation region are made in Sec. VI. The conclusion is given in Sec. VII. 



II. VALON AND PARTON DISTRIBUTION FUNCTIONS 

The distribution functions of valons and partons have been considered in Refs. They were, however, 

determined by fitting the muon and neutrino scattering data of the late 70s. We now have modern parton distributions 
from various groups at various values of . It is therefore appropriate for us to revisit the problem and determine 
the valon distribution functions in light of the new parton distribution functions. 

In the valon model a proton is considered to consist of three valons {UUD), which have the same flavors as the 
valence quarks [uud) that they individually contain. Thus a valon may be regarded as a parton cluster whose structure 
can be probed at high Q^, but the structure of a nucleon itself in a low-p-p scattering problem is described in terms of 
the valons. As in the parton model we work in a high- momentum frame so that it is sensible to use the momentum 
fractions of the constituents. Reserving x for the momentum fraction of a quark, we use y to denote the momentum 
fraction of a valon. In this paper y never denotes rapidity. Let the exclusive valon distribution function be 

GuuD{yi,y2,y3) = 9 {ymTys S{yi + 2/2 + y3 - 1), (i) 

where yi and y2 refer to the U valons and y^ the D valon. The delta function ensures that the three valons exhaust 
the momentum of the proton. The exponents a and (3 will be determined by the parton distribution functions. The 
normalization factor g is determined by requiring that the probability of finding these three valons in a proton be 
one, i. e. 

dyi I dy2 dy3GuuD{yi,y2,y3) (2) 



Thus we have 

g = [B{a + 1,(3 + l)B{a + I,a + f3 + 2)]"^ , (3) 
where B(m,n) is the beta function. The single valon distributions are obtained by appropriate integrations 

Gu{y)= ! dy2 I dy3GuuD{y,y2,y3)^gB{a + l,P+l)y''{l-yr+^+\ (4) 



Goiy)^ J dyi J dy2GuuD{yi,y2,y)^9B{a + l,a + l)y'\l-yf'-+\ (5) 

The two- valon distributions are trivial because of the (5- function in Eq. (|^). 

In a deep inelastic scattering the structure of the proton is probed to reveal the parton distributions, which in the 
valon model are convolutions of Gy{y) with the evolution functions that describe the valon structure. The latter have 
two varieties: K{z, Q^) for the favored partons and L{z, Q^) for the unfavored partons. That is, u is favored in U and 
d is favored in D, but they are unfavored in D and U, respectively. At high we may ignore the influence of the 
spectator valons when one valon is probed (the usual impulse approximation), so we can write for the u and d quark 
distribution functions as 

xu{x,Q^)^ [ dy[2Gu{y)K{x/y,Q^) + GD{y)L{x/y,Q^)], (6) 



xdix,Q^)^ I dy[GD{y)Kix/y,Q^) + 2Guiy)Lix/y,Q^)] 



(7) 
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It should be emphasized that u{x), d{x) and G{y) are noninvariant distributions defined in the phase space dx and 
dy, while K{z) and L{z) are invariant distributions defined in the phase space dz/z. 

The favored distribution K(z,Q^) has two parts, valence and sea, while the unfavored distribution L{z,Q'^) has 
only sea. The valence part is also referred to as the non-singlet component, Kns{z, Q^), so we may write 

K{z,Q^)^Kns{z,Q^)+L{z,Q^). (8) 

At high both K^s and L can separately be determined in perturbative QCD. However, since low-Q^ parton 
distributions are now available by extrapolation, and especially since we shall apply the valon model to low-px 
processes (and therefore at low Q^), we use phenomenological forms for K^s and L at low with parameters to be 
determined by the low-Q^ parton distributions. We adopt the forms 

KNs{z) = z''{l-z)''/B{a,b+l), (9) 



L(z) =4(1-^)', 

where Eq. satisfies the requirement that there is only one valence quark in a valon, i. e. 

f'^K^siz)^!- (11) 
Jo z 

Eq. ( |l0| ) has the usual sea-quark distribution, also used previously [l^ . 

The parton distributions that we use to fit are the ones determined by the CTEQ collaboration [ p^ . In particular, 
they have the distributions at low Q^, labeled CTEQ4LQ, posted on the web We choose the one a,t — 1 GeV^ 
evolved from Qq = 0.49 GeV'^. We fit the u{x, Q"^) and d{x, Q"^) distribution functions using Eqs. (P)-(p^) by varying 
a, /?, a, h and l^- The results are shown in Fig. 1, where the CTEQ functions are in solid lines, and our fitted curves 
arc in dotted lines. We have attempted to fit the u-quark distribution as perfectly as possible; the fit of the d-quark 
distribution turns out to be good only at high x. Missing the normalization of d{x, Q^) at low a; is a blemish, but is 
acceptable since we shall use the valence quark distributions mainly in the large-x region. Besides, the reliability of 
an extrapolation of high-Q^ deep inelastic scattering data to low can always be called into question. 

The parameters of the fit are 

a = 0.70, /3 = 0.25, (12) 



(10) 



a = 0.79, 6 =-0.26, 4 = 0.083. (13) 

The values of a and 13 are not very different from the ones determined in Ref. pl based on muon-scattering data 
at — 22.5GeV^ analyzed by Duke and Roberts ||l8|; there we had a = 0.65 and P = 0.35. As can be seen from 
Eqs. (||) and (||), our present result means that Gu(y) (1 — y)^'^^ and Gniy) oc (1 — y)^-^ at large y, as opposed 
to (1 — y)^'^ and (1 — y)^'^, respectively, for the previous result. The values of a and b in Eq. (^3|) suggest that 
KNsiZjQ"^) is highly peaked near z — 1, according to Eq. (^). That is as it should be, since is only about twice 
the value of Qq and thus not much evolution. If there were no evolution, i. e., = Qq, then Kns(z,Qq) would be 
S{z — 1) and L{z) = 0, a situation indicating no probing of the valons. What we have at = 1 GeV^ gives only a 
modest degree of resolution of the valon structure, and K]s[s{z,Qq) is changed from 5(z — 1) to (1 — z)^°'^^ in the 
large y region, while L{zj Qq) becomes non- vanishing. These are the valence and sea quark distributions that we shall 
use for nucleon production at low pr in the following. 



III. MOMENTUM DEGRADATION IN THE VALON MODEL 



The valon model for inclusive reactions is basically a s-channel description of particle production in contradistinction 
from the Regge-Gribov approach jl^ which is essentially based on cutting t-channel exchanges in elastic amplitudes. 
The two approaches are roughly complementary in that the former can best describe the fragmentation region, while 
the latter is more suitable for the central region. 

It has been known for a long time that the pion inclusive cross section in the proton fragmentation region has a 
xf dependence that is very similar to the quark distribution in the proton pof . That similarity has been reconfirmed 
more recently in pA collisions by the E910 experiment even at the relatively low energy of E^^^-^ = 17.5GeV. It 
suggests that the proton structure is highly relevant to the spectra of particles produced in the fragmentation region. 
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The connection between the quark distribution and the pion inclusive distribution was first put on a calculable basis 
by the recombination model jl^ , which was subsequently improved in the framework of the valon model Jl2[ . The 
idea of recombination as a basis for hadronization in hadronic collisions is eminently reasonable, since for a pion to 
be detected at xf — 0.8, say, it is far less costly to have a quark at xi coalescing with an antiquark at X2, each < 0.8, 
to form a pion at x^? = xi + a;2, as compared to the fragmentation process for which a quark or diquark must first 
have X > 0.8. In quark jets, fragmentation is reasonable because hadronization is initiated by a leading quark, but it 
is less persuasive when adapted to hadronic fragmentation unless the projectile momentum resides entirely in a quark 
or diquark j2l], contrary to Ochs' observation of the relevance of the canonical quark distributions. 

The s-channel treatment of relating quark distributions to inclusive hadron distributions in pp collisions essentially 
regards the effect of the opposite-going initial proton as unimportant, an approximation that can only be justified in 
the fragmentation regions due to short-range correlation in rapidity. Thus the valon model, as it has been developed 
up to now, is not expected to be applicable to the central regions, where the interaction between the two incident 
particles is of paramount importance. For pA collisions even the factorization of the fragmentation properties is not 
entirely valid, since the hadron distributions are known to depend on centrality or target size. This is the problem 
that we shall treat in this paper in the framework of the valon model. More specifically, we shall consider the problem 
of momentum degradation of the valons due to the nuclear medium. Recombination occurs outside the nucleus and 
is therefore unaffected. 

The NA49 data on p-Pb collisions are presented in terms of two mean values, D, of the number of participating 
nucleons in the target: D = 6.3 for central collisions and 9 = 3.1 for non-central collisions We assume a Poissonian 
fluctuation from that mean with the distribution 

P,(^) = !^(e^-l)-i, (14) 



which is normalized by 



(15) 



where we have excluded the ;/ = term, since it is necessary for v > 1 in order to have a collision. Thus v is the 
number of nucleons in the nucleus that suffer inelastic collisions in any given event, v being an integer. That is the 
counting on the target side, while on the projectile side we count in terms of the valons. Let the ith valon encounter 
i/i collisions. In general, the total valonic collisions is bounded by 

3 

<'^iyi <3v. (16) 

i=l 

The upper bound occurs only when all three valons participate in each of the struck nucleon, while the lower limit is 
for only one valon per struck nucleon. In the Appendix we shall show that the data favor the lower bounds, so for 
simplicity we proceed in the following with the assumption 

v = + 1^2 + vs- (17) 

It is useful to interpret this in the i-channel picture. The incident proton is represented by three constituents, each 
of which exchanges Vi ladders with Vi target nucleons so that the overall diagram is highly non-planar. Cutting the 
ladders gives rise to the particles produced in the s channel, mostly in the central region. 

Focusing on the evolution of the valons, we first assume that the three valons interact with the target independently, 
since they are loosely bound to form the proton, just as the nucleons are in a deuteron. If we denote the degradation 
effect of the nucleus on the ith valon by D{zi,iyi), then we can write the evolution equation on the valon distribution 
as 

y'iyW3G'iy'i,y'2,y3;'^i,''2,>^3)^ f dyidy2dy^G{y^,y2,y3)D{y^,vi)D{y^,v2)D{y^,v^) (18) 

J yi y2 ys 

where G' is the valon distribution function after (i^i, 1^2, t's) interactions with the nucleus. As in Eqs. (^ and (|^), 
G and G' are noninvariant distributions defined in the phase space dyi dy2 dy^, etc., while D{zi,Vj) is an invariant 
distribution defined in dzij Zi. We postpone our discussion of what D(zi, Vi) is until Sec. V. 

For an event with v collisions the evolved valon distribution function G'j,(2/^, 2/2; 2/3) with v partitioned as in Eq. 
( p^ ) is given by the multinomial formula 

G'Ay'i,y'2^y'3) = ^E \^ \ . G\y[,y'2,y'3;vi,iy2,i^3), (19) 
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where [ui] implies that the summation is over vi^V2, and subject to the constraint of Eq. (17). Note that if the 
nucleus had no effect on the valons, i.e., 

D{z„v,)=5{z,-\) (20) 

independent of Vi, then G'{y'^, y'2, y'3; vi,V2-, 1^3) becomes G{yi,y2,yz) in Eq. (|l|) and so does G"^(yi , ^3) ^ Eq. (|l|), 
as it should. To relate to the experimental i>, we have 

G'M,y2. 2/3) = E G',{y[,y'2,y'3)P,{u), (21) 

where Piy{i^) is given in Eq. (|T^). 

It is useful to introduce the notion of an effective nucleon after v collisions by giving it a momentum fraction y' 
and defining the probability of finding it at y' by 

QUv') = I dy[dy'2dy', G'M,y'2,y',) Siy[ + + - y') (22) 
Conservation of baryon number requires that 

dy'Gliy') = 1=1 dy[dy',dy',G'M,y'2,y',). (23) 



The possibility that flavor can change is a secondary issue that will be discussed later; here, the issue is that the 
baryon should be somewhere in the interval < y' < 1. The second half of Eq. (|2^) puts a constraint on D{zi,i'i), 
since we can obtain from Eqs. 1^ and (|) 

-Diz,u,) = l. (24) 

z 

The likelihood that the three evolved valons will in reality reconstitute a nucleon is extremely low, but the fictitious 
nucleon that they form carries a baryon number that is conserved, and a momentum that is not conserved. Indeed, 
we expect the average y' to decrease with i^, i. e., 

y'.^ C dy'y'gl{y')<l. (25) 

That is commonly referred to as stopping. In Sec. V we shall infer from the data what the stopping power is. 

Even without stopping, such as in pp collisions, it does not mean that the real proton produced cannot have xp < ^■ 
It is known that in pp collisions the proton inclusive cross section da/dxp is nearly flat in xp- Stopping goes on top 
of that distribution, making it roughly exponential decrease in xp- How to proceed from the valon distribution G' to 
the detected proton distribution Hp is the subject of the next section. 

The convolution equation (^) can be simplified when expressed in terms of the moments on account of the convo- 
lution theorem. Thus let us define 

/■I 

D{ni,Vi)^ / — -z'l'^^D{zi,Vi), (26) 



G(ni,7i2,"3) ^ dy 



dy2 



1-2/1-2/2 



dy3 



n 



Ui-l 



G{yi,y2,y3) 



and similarly for G'{ni,n2,n3; 1/1,1/2, 1^3) and G'^{ni,n2,n3). It then follows from Eqs. (^8|) and (p^ ) that 

1 iy\ ^ 
G'U"i,«2,«3) = ^y]~i ''\ :G{ni,n2,n3)Y\ D{n^,iy,). 



(27) 



(28) 



Now, D{ni, Vi) itself can be described by a convolution equation If instead of the discrete Vi we use a continuous 
variable L that denotes the length of the nuclear medium a valon traverses, we can express the change on D{z, L) for 
an incremental distance dL in the form E3| 



d 

dL' 



dz' 

D{z,L)^ / —D{z',L)Q{z/z') 

J z Z 



(29) 
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with some reasonable kernel Q{z/ z'). In terms of the moments with 

Q(n)= t dCC'^QiO. 
Jo 



Eq. ( p^ ) becomes 



whose solution is 



The constraint (EJ) implies 



-^D{n,L)^ Din,L)Q{n), 



D{n,L) ~ exp Q{n)L 



D{1,L)^1, and Q(l) = 0. 
Since Vi is proportional to L, let us now revert D{n, L) to D{ni, ui) and write Eq. ( |3^ ) as 



(30) 



(31) 



(32) 



(33) 



(34) 



where d(ni) is trivially related to e'^^"'^ with a power exponent whose detail need not be specified here. From Eq. 
(P3|) follows 



d{l) = 1. 

We now can use Eq. (jsj) in ( p8| ) and obtain 

G"^,(ni,n2,n3) = G(ni, 712,713) 



1 ^ 



(35) 



(36) 



What we have derived here is that the dependence on v is in the exponent, implying that the effects of the successive 
collisions with the nucleons in the target nucleus are multiplicative, as is reasonable. We can now go back to Eq. ( [25| ) 
and calculate the average yl after collisions. Using Eq. (|2^) we get 



yl = GU2,1,1) + GU1,2,1) + G;(1,1,2) 

= G(2, 1, 1) + G(1, 2, 1) + G(l, 1, 2)] {[2 + d(2)]/3r 

with the help of Eqs. (BSh and (pq). The first factor in the square brackets is just 1, since it is 



which is guaranteed by the (5-function in Eq. Hence, we have 

y:, = r, e-[2 + d(2)]/3<l. 

Averaging over Py[h'), defined in Eq. (jl^), yields 

00 



(37) 



(38) 



(39) 



(40) 



This is the average momentum fraction of the effective nucleon after v collisions but before fragmentation into final- 
state particles in the fragmentation region. 
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IV. FRAGMENTATION AND RECOMBINATION 



We now consider the problem of how a projectile proton fragments and how the quarks recombine to form the 
detected nucleon, thereby specifying how the inclusive distribution can be calculated. To give an overview of the 
procedure, let us summarize the two steps above by the two invariant distributions: X2, 0:3), the probability of 

finding a u quark at xi, another u quark at X2, and a d quark at X3, and Rp{xi,X2,x^, x), the recombination function, 
which specifies the probability that those three quarks coalesce to form a proton at x. How these distributions are 
related to the valon distributions will be discussed later. But first we state that the invariant distribution function 
for the detection of a proton at x is 



X da^ 1 /" dxi dx2 dx^ 
r~ = Hpix) = — I 

(Tin dX N J Xi X2 X3 



^p(^) = Tr I F{xi,X2,X3)Rp{xi,X2,X3,x) (41) 



where the normalization factor N will be given below. Eq. (^) is the essence of the recombination model ||T^-||lJ|. For 
meson production, only the distributions for q and q need be considered and Eq. ( pl| ) can be simplified accordingly. 
As in all distributions considered in this paper, color and spin components are averaged over in the initial state and 
summed in the final state so that we are not concerned explicitly with such degrees of freedom. Flavor, however, is 
different, since we identify the final-state particles by their flavors; that problem will be treated presently. 



Before proceeding we emphasize what has been mentioned in the preceding section already, namely: Eq. (41) is 
expected to be valid in the proton fragmentation region only, if F(xi,X2,X3) is to be determined from the projectile 
valon distributions with no quarks originating from the target nucleons. In this s-channel approach, the factorization 
of the projectile and target fragmentations, apart from the momentum-degradation effect studied in the previous 
section, can be justified only if the two fragmentation regions are well separated. At AGS that is not the case. Even 
at SPS the central region in rapidity can encompass sizable portions of the positive and negative xp variables. The 
application of the valon model to the analysis of the data therefore needs some help from the experiments. 

Fortunately, the NA49 collaboration has treated their data in such a way as to eliminate the contribution of 
the target fragmentation from the projectile fragmentation region. From their data on the net proton produced, 
p + A ^ (p — p) + X, for which we use the abbreviated notation (p — p)p, they subtract the distribution for [p — p)t^, 
which is ^[{p — p)tt+ + {p — p)7T~]- By charge conjugation symmetry, {p — p)-^- should have no projectile fragmentation, 
only target fragmentation. Thus the difference {p — p)p — {p — p)tt should have no target fragmentation |2^. With 
those data as our goal for analysis, Eq. (^ ) is then particularly suitable. 



The normalization factor N in Eq. (41|) is determined by requiring that the proton distribution in the absence of 



es 



/''-770(x)= [\xH°p{x) = l, (42) 
^ Jo 



any target, Hp{x), satisfies the sum ru 



which follow from the condition that the number of proton and its momentum fraction be 1 . Without any collision 
the quarks are identified with the unresolved valons, so F{xi,X2,X3) becomes 

,X2,X3) = XiX2X3GuUd{xi,X2,X3). (43) 

The recombination function is the time-reversed form of the valon distribution, so 

R ^ ^ ^1^2X3 ^ (Xi X2 X3\ 

Rp(Xi,X2,X3,X) = ^ — GuuD — , — , — ■ (44) 

x-^ \ X X X J 

Putting these in Eq. (p]), we obtain in view of (|l|) 

jyO(x) = C /' dx^ C dx2{xrX2f^+\l ~ x^ - X2f^+^x-'^^^+^+^^{x - 1). (45) 



N Jo JO 

Because of the presence of 5{x — 1), the two integrals in Eq. (^) are identical, and we get 

N = g^B{2a + 2,2a + 2/3 + 4)B(2a + 2,2/3 + 2). (46) 

The factor g^, although known from Eq. (|^), will cancel the similar factor that will emerge from the integral in the 
numerator of Eq. (p]), just as they appear exphcitly in Eq. (^. 

The identification of the recombination function with the invariant form of the valon distribution in Eq. (^4|) is the 
principle characteristic of the valon model. On the one hand, it recognizes the role of the wave function of the proton 
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both in a projectile and in a produced proton. On the other hand, the momentum fractions Xi of the outgoing valons 
can add up to a proton at x, so there is no need for any constituent in the process to have a momentum fraction 
greater than x, as would be necessary in a quark fragmentation model. One may then question how in a collision 
process can the quarks at xi,X2, and X3 in Eq. ( p] ) become the valons of the outgoing proton. The answer is that 
hadronization occurs outside the target, and that the quarks moving downstream dress up themselves and become 
the valons of the produced proton without any change in the net momentum of each quark/valon, which is all that 
matters in the specification of Rp{xi,X2,X3,x). 

We now consider the quark distribution i^(a;i, a;2, a::3) in Eq. ( |4l| ) in pA collision before recombination. In the 
preceding section we have formulated the procedure to calculate the effect of the nuclear medium on the momenta of 
the valons as they traverse the target. Momentum degradation is, however, only one of the effects of valon-nucleon 
interaction. If such an interaction is represented by a Regge exchange, we should also consider the possibility of 
flavor changes of the valons due to non-vacuum exchanges at non-asymptotic energies. In the spirit of the s-channel 
approach that we have taken, in which the probabilities of occurrences at various stages are assembled multiplicatively, 
we assume that the flavor changes at each of the I'i collisions are incoherent so that the net probability of a flavor 
change after i^i interactions is also multiplicative. Let q be the probability of a flavor change from U to D, or from 
D to U, at one of the Vi interactions. Furthermore, let q^. be the probability of flavor change after Vi interactions. 
Then q^. satisfies the recursion relation 

Qi^i+i ^ q-y,{^ - q) + - q-y,)q, (47) 

where the first term on the RHS denotes no change in the last step from to Vi+i, while the second term denotes a 
change in the last step. The solution of Eq. ( ^7| ) is 

q..=l[l-{l-2qr]. (48) 
We may now write what a U and a D valon become after vi interactions in obvious notation: 

U^p,^U + q,^D (49) 

D^p,^D + q,^U (50) 

where p^. = 1 — g^. . This regeneration process depends on one parameter q. We expect q to decrease with energy ^/s. 
Here we treat it as one free parameter to fit the NA49 data at one energy E-^^-^ = 158 GeV. 

For the quark distribution in a valon we have the favored and unfavored types discussed in Sec. II, and denoted 
by K{z) and L{z), respectively. We drop the dependence, since we now consider \ow-px hadronic processes for 
which there is no precise Q^. Nevertheless, we shall use the parameterization in Eqs. (^, ( p^ and (p^), from the 
— 1 GeV^ CTEQ parton distributions. We use the following notation for the invariant distributions of quarks in 
valons with superscript / signifying 'favored' and u 'unfavored': 

V/j: favored quark at xj in valon at y'^, 

wfjf.: two favored quarks at xj and Xk in valon at y'^. 

Other distributions involving unfavored quarks are similarly defined. Examples of V/j, Fj", W/jf., and W^/"^. are 
depicted by diagrams in Fig. 2. 

In view of Eq. (p9[) and Eq. (&□) we can write by definition 



(51) 



(52) 



/ )jk^ ^''^ [\y't)^ \y'i - ^3 ) )jk ^^^^ 



where 



/i ( 3 ) 



Vi 



Xk 



(56) 



In terms of these V and W distributions we can now write out, by inspection, aU possible contributions to the uud 
quarks, shown in Fig. 3, for the production of a proton 



+^yi,w^h + 2{V3"ii^i%}i2 + '^viwii^, 



(57) 



where { }i2 denotes symmetrization of xi and X2. We have ignored the contributions corresponding to all three quarks 
coming from the same valon. The quark distribution from the proton source for proton production is then 



(58) 



Fp{xi,X2,xs,) = I dy^dy2dy^Gp{y-^,y2,y:i)Mp{y^,y2,y^;xi,X2,X3) 



where G'i^{y[, 2/2: J/a) given in Eq. (pl|). For p production we need only change all K functions in Eqs. (pl|) - ( pq ) to 
L functions, since all anti-quarks are in the sea. Denoting the corresponding quantities in Eqs. (|57| ) and (p8|) by Mp 
and Fp, we have finally for net proton production 



Fp.p = ® [Mp - Mp), 



(59) 



where the convolution is defined by the integral in Eq. (|58|). It should be recognized that Gp(?;^ , ?/2j 2/3) involves a 
summation of G'{y'i, y'^) over which in turn involves a summation G'{y[, j/jj ^3! '^Ij ^2, 1^3) over vi, U2 and uy, that 
appear in the V and W distributions. 

We are now ready to substitute F(xi,X2,X3) and R{xi,X2,X3,x), defined in Eq. (^), into Eq. ( |4l| ) to calculate 
H{x). Nine convolution integrals are involved: yi,y[, and Xi. Obviously, we should go to the moments and reduce 
them to products. First, using Eqs. (0) in (44), we have 



H{x) = Aa;-(2o+/3+2) / dxidx2dx^F{xi,X2,X^){xiX2T+^X^^+^ 5{xi + X2 + x^ ^ x). 



For convenience, let us leave out the known factors and define 



H'{x) 



N 



-X 



2a+/3+4 



Hix) 



(60) 



(61) 



Then define the moments 



H'{n) = / dxx'^-'^H'{x), 
Jo 



(62) 



Thus it follows 



F{ni,n2,n3) = J dxidx2dx3 ^^^^x"' 



F{xi,X2,X3). 



ni\n2\n'i\ 



F[ni + a + 3, 712 + a + 3, 713 + /? + 3). 



(63) 



(64) 



From Eq. ([Sq ) we expect from the convolution theorem to have 

F(ni, 712,713) ^ G"p(7ii,ri2,ri3)M(7ii, 712,773) 



(65) 
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where G"p is given by Eq. (p8[), apart from the Poissonian sum of (^). However, because of the ordering of xj relative 
to y'^ in Eq. (p?]), the simple form of ( |6^ ) is valid only for the first term of M{y[, j/j? Ua': 2^3)- For that first term, 

which we denote by F'^^\ we have 



F^'^\ni,n2,n^) = G"(ni, n2, ns) 



(66) 



where the subscript u has been omitted. The summation over Vi in Eq. (^8|) for G'(ni, 712, na) should extend over 
and q^. in Eq. (pq). The moments of K{zi) and L{zi) are defined as usual for the invariant distributions as 



(67) 

For all other terms in Eq. (^7| ) there are and functions, and the simple form of Eq. (^5|) must be modified, 
especially when only two valons contribute. For notational simplicity let us define 



ff 



i,jk 



- p^.{K{nj,nk)L{nk)}jk + quA^inj ,nk)L{nk)} jk, 
■ p^.{L{nj,nk)L{nk)}jk + q,y^{K{nj,nk)L{nk)}jk, 



where K{nj,nk) is defined by 



(68) 
(69) 
(70) 
(71) 
(72) 

(73) 



and similarly for L{nj, Uk). These moments arise whenever two quarks are from the same valon, as they do for all W 
distributions. Take, for example, the third term in Eq. (p7h; we have 



dv'2 G'jju{v'i,v'2) 



' Vi ¥2 ¥2-^2 



dy[ 



dy'2G'uuiy'i,y2)yT'-'yT'^'''-'2p^^p^, 



dziz'^'-'Kizi) / dz2Z^'-'{l^ Z2r'-'Kiz2) / dZ3Z^'-'L{z3) + 
Jo Jo 

= 2{G'jjjj{n,,n2 + - l)l//iW^2%}i2, 
where {}i2 here means symmetrization of ni and 712, and 

G'uui^li'^2) = G"[j[/£, (mi, 7712,7713 = 1). 

Performing the same type of operations on all terms we obtain 

Fp(77l,772,773) 

= G'uuD{ni,n2,n3Mvlvl + 2{G'uuDini,n3, n2)V/,V,"2} 12V2I 



(74) 
(75) 



+ 2{G'^^(77l, 772 + "3 - l)VAWi;i,},2 + 2G'^^(773, 77l + 772 - l)V,lW^;i2 



+ 2{G'uDini,n2 + 773 - l)^AW^3,32}l2 + 2G'c,^(773 , 77i + 772 - l)Fi"3W^3^^2 

2{G'uD{n2 + - 1, «i)W^/23^3"i}i2 + ^{G'uDini +n2-l, r73)M^i%V^33}i2- 



(76) 
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Substituting this in Eq. (p4|), we have the final form for the moments Hp{n). For p production we need only change 

all K in Eqs. (||) - (H) to L. 

For the production of neutron the recombination function is 

Rn(Xi,X2,X3,x) = — Gddu — , — , — , (77) 

X'^ \ X X X / 

where we use j/i and 1/2 to refer to D and j/3 to U, so the dependence of 0001/(1/1,1/2,1/3) on yi is the same as that 
of G(7[/£)(yi, 2/2, ys) given in Eq. (|l|). The quark distribution Fn{xi,X2,X3) with the corresponding identification of 
xi and X2 with the d quark, and X3 with the u quark, also has eight terms, as shown in Fig. 4. We can write, by 
inspection, the moments Fn{xi,X2,X3) similar to Eq. ( [76| ) 

F\{Ci,X2,X3) = 

HG'uuDin3,ni,n2)V/,V,\V^,h2 + G'^c,c(ni, ^2, ns)V,\V^-,V^l 
+2G'uuin3, n,+n2~ 1)^/3 W^2"i2 + 2{G'uuini,n2 +713- l)n'!W^2^32}i2 
+2G'uoin3, ni + ^2 - 1)^/31^3^(2 + 2{G'y^(ni, 712 + ^3 - 1) Vi"iW^3%}i2 

+2G'uD{ni +n2~l, n3)Wi';,2V3'',3 + HG'uoini + "3 - 1, «2)W^/3iV^3-^2}i2- (78) 

Substituting this in Eq. (|6^ ) gives H'^{n). For n production simply replace all K in Hn{n) by L. 
Finally for net nucleon production we have 

ij;_p(n) - H'^in) - %{n), (79) 

K^^{n)^HUn)~Hk{n). (80) 
We shall use these in the next section to determine the x dependences to be compared with the data. 

V. NET PROTON AND NEUTRON DISTRIBUTIONS 

Before we can compute the distributions H{x), we need to specify quantitatively the degradation function D{zi, Vi) 
introduced in Eq. ( p^ ) and discussed between Eq. ( |29| ) and (|35|). We proposed an evolution equation for D{z,L) in 



Eq. (29) but left the kernel Q{z/ z') unspecified. Now, to proceed we must specify (5(C), which is uncalculable because 
it represents the non-perturbative effect of the nuclear medium on a valon as it propagates an incremental distance. 
We shall use an one-parameter description of the effect, so we shall be approximate by assuming that the effect is like 



a one-gluon exchange ||l^, i. e. 



(1-0 

where the singularity at C = 1 is regularized by the subtraction 



QiO = 7T^^, (81) 



1 1 . f'^ dx 



^(C-1) ^ • (82) 



(1-0+ 1-C 'J^ i-x 

Evidently, Eq. ( |8l| ) satisfies the condition 

' ^Q(0 = 0, (83) 
s, 

which is required by the constraint Q(l) — stated in Eq. (^Sh that follows from baryon conservation. From the 
definition of the moments given in Eq. (pG) we find, using Eq. (|8l|), 

n-l 

Q(n) = -At^ = -«[^(n) + 7£], (84) 
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where ipin-) is the digamma function and 7^ is the Eulcr's constant, 0.5722. Substituting Q{n) into Eq. (|3^), and 
changing L to Vi that involves a constant factor, thereby effecting a change from one unknown parameter k to another, 
k, we obtain the form for D{ni, v^) in Eq. ([sj) with 

d{ni) ^ exp{-fc[V'(ni) + 75]}. (85) 

This is a one-parameter description of the effect of momentum degradation. We shall vary k to fit the data. Eq. ( |85| ) 
is a rigorous consequence of a simple form for Q(C) given in Eq. (|8l|), whose reliability is unknown. The validity of 
d{ni) as expressed in Eq. ( ^5| ) can only be inferred a posteriori from the fit of the data. The exponential dependence 
on the degradation strength k follows only from the linear dependence of Q{C) on k, and is sensible. 

We have only two free parameters, k and q, to vary to fit the NA49 data on p-p and n-n ||] . Recall that q is introduced 
in Eq. ( ^ ) in connection with flavor changes. We repeat that the data do not include target fragmentation because 
{p — p)-n and {n — n)^ have been subtracted out. Thus the data represent only proton fragmentation and are ideal for 
our analysis by Hp-p{x) and Hn-n(x). 

Since the data are in the x variable, we must make the inverse transformation from our moments to H{x). Instead 
of making the inverse Mellin transform, let us exploit the orthogonality of the Legendre polynomials and shift the 
variable to the interval < x < 1. Thus, define 

giix) ^ Pe{2x - 1) (86) 

so that 

1 1 

dxgi{x)gm{x) = 2f + l *^^"'' ^^"^^ 

If we expand the distribution H'{x) in terms of gi{x) 

00 

H'{x)^Y.^2l+\)htgt{x), (88) 
£=0 

then the inverse is 

hi^ f dxH'{x)ge{x). (89) 
Jo 

These hg can be expressed in terms of the moments H'{n) if we express gg{x) as a power series in x 

i 

5Ka:) = ^a^a:^ (90) 
where a\ are known from the properties of Pi{z). Thus from Eq. (p9) we have 



hg = ^a^pH'{i + 2), (91) 

1=0 



where H'{n) is defined in Eq. (|62|). It is now clear that our theoretical results in H'{n) can be transformed to H'{x) 
through Eqs. (^8|) and (^) once we have the coefficients a\. Furthermore, if H'{n) becomes unimportant for n > N, 
then the sum in Eq. (p8|) can terminate at N. 

To determine a^, we make use of the recursion formula 

(^+ l)P,+i(z) = (2£+ l)zP,(z) -£P,_i(z) (92) 

to infer through Eqs. ( ^6|) and ( |90| ) 

a? = -i[(2£ ~ l)aU + {£- l)aU], (93) 

4 = -j[{2i - l){aU - '2a]r\) + - 1)«L2], (94) 
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where £ > 2, and 1 < i < i. For i — and/or £ < 2, we have 

a[J = 1, a? = -1, a\=2. (95) 

With these we can generate all a], so hi can be directly computed. 

Summarizing our procedure, we calculate Hp_p{n) and Hj^_f^(n), substitute them in Eq. (^Tj) and then (^8|), and 
then use Eq. ( pi] ) to determine Hp-p{x) and Hn-n{x)- We vary k and q to fit the data of NA49 shown in Fig. 5. The 
inclusive distribution dN/dxp corresponds to our H{x)/x. The solid lines are our results for p ~ p and the dashed 
line n — n for both v — i.l and 6.3. The values of the parameters adjusted are 

k = 0.62, q = 0.37. (96) 

The most striking aspect of our result is that the normalization of the calculated distributions turns out to be correct, 
even though we have no free parameter to adjust that. The degradation strength k affects the shape of the distributions 
and the flavor-flip probability q affects the difference between p — p and n — fi. The agreement between theory and 
experiment is fairly good, considering that we have only two free parameters and that the experimental errors (a 
typical size of which is shown in the figure) are large, especially at large x. The shapes of the distributions for p ~ p 
are reasonably well reproduced and so is the dependence on D. For n — h the calculated curves are somewhat steeper 
than the data. However, there exist some data points for n — n above x = 0.8 that are significantly lower, though with 
much larger errors. Taken as a whole the agreement is satisfactory. Thus we conclude that the physical process of 
proton fragmentation and the nucleon momentum spectra are well understood in the framework of the valon model. 
We can determine from the value of k what the momentum degradation length is. First, we get from Eq. ( |85| ) 

d{2) ^ e-'' = 0.54 (97) 

for k = 0.62, since '4>{2) + 7^; = 1. By definition Eq. (j3^) we obtain ^ — 0.85. Inserting this in Eq. ( pO[ ) we obtain a 
dependence of {y')p on 9 that can be well approximated by 

(2/% cx e-(l-^)^ (98) 
for V >2. This gives the fractional momentum loss per collision 

If we related v to average nuclear path length Lhy u — apptA and Ia = ApL, where p = ^i?^, then for app = 30mb 
we have D ~ 0.4L with L in fm. Thus, if we define the degradation length A by 

{y') = e-^/^, (100) 

then 

A = [0.4(1 -0]"^ - 17fm. (101) 

This gives an estimate of how far a proton must travel in a nuclear medium in order to lose its momentum by a factor 
of e-i. 

The value of g = 0.37 for flavor-flip probability may at first sight appear to be surprisingly large. However, if it 
is regarded as an effective way of accounting for resonance production, it becomes quite acceptable. To see that, we 
first state that resonance production, which we have not taken into account explicitly, can easily produce neutron 
from \uud) through A"*" n + Tr~^ . The process can be depicted by a dual diagram, as shown in Fig. 6(a). Such a 
leakage of -I- charge through the emission of 7r+ is equivalent to a flavor flip, which changes UUD to DUD, symbolized 
by a square box in Fig. 6(b), and the favored process, for example, of having valence quarks changes from uud to 
dud. Since a consideration of resonance production would involve masses, threshold, polarization, decay distribution, 
and other complications, our method of using flavor changes to account for the effect presents considerable technical 
economy. Even though resonance production can occur only at hadronization in the end charge leakage can take place 
at any point where a projectile valon interacts with the target; hence, the consideration in Sec. IV leading up to Eq. 
( pO| ) is an effective way to take such subprocesses into account. 
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VI. PION PRODUCTION 



Having successfully computed the nucleoli distributions, resulting in the determination of the only two free param- 
eters in the model, we are now able to predict the pion distributions without any further ambiguities. The NA49 
data presented in Ref. Q do not include the pion spectra. The E910 data [D do have the pion distributions in the 
proton fragmentation region; however, being at E^^-^ = 12 GeV the energy is too low to avoid substantial spill-over 
of quarks and produced hadrons from target fragmentation into the a; > region. Without the target fragmentation 
being subtracted, as is done for p — p in the NA49 data, the E910 data cannot be compared to the predictions from 
our model. We present our result below for a future comparison. 

The invariant distribution for pion production is 

- H^{x) = / — — F^(a;i,a;2)i?7r(a;i,a;2,a;) (102) 



nn 



dx J X\ X2 



where Ft^(x\^x-2) is the invariant distribution for finding a quark at x\ and an antiquark at X2^ and -R^r is the 
corresponding recombination function to form a pion. An important aspect about pion production concerns the role 
of gluons, a subject that will be discussed below in connection with L(z). For now, we consider the appropriate forms 
for and Rt^. 

We begin with the valon distribution G'i^{y[,y2, y'^), as given by Eq. (^l|). As before, we shall omit the subscript P, 
and replace it by the valon labels so that for the two-valon distributions we have 

G'uu{y'i.V2)-^ I dy',G'uuD{y'i,y2,y'3), (103) 



and similarly for G'jjjj{y'i,y'^) by integrating out j/j- The single- valon distributions are involved in the 5 subprocesses 
that contribute to F.^{xi, X2) shown in Fig. 7 for the production of 7r+. Of course, because of the flavor changes, the 
valon labels are only indicative of the unchanged components. More precisely, we can express F^^ in the moment form, 
as in Eq. ([76|), for 7r+ 



F,+ (ni,n2) = [G'yy(ni,n2) + 2G'^^(ni,n2)]{F/iF2"2}i2 
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+2G'uD{n2,ni){V,'i,V2\}i2 + 2G'jj{ni +n2- l)Wl 
+2Gi3(ni+n2-l)t^3^^2, (104) 



and for tt 



F,-(ni,n2) = G'uu{'^i.ri2)vlV22 

+2G'j;^(n2,ni){f/iFi"2}i2 + 2G'jjj,{nun2){Vl^V^2}i2 

+2G'jj{m +n2- l)W^^i2 + '^G'oim +n2- l)W'/"2- (105) 
For the recombination function we have 

U ( \ ^1^2^^ (Xi X2\ , , 

Rt,{xi,X2,x) ^ —Gjjq\—, — \, (106) 

X ^ X X ' 

where G'^jj{yi,y2) is the valon distribution in a pion. For the latter we adopt the same form derived in Ref. 

GlBiyi^y2) = S{yi+y2-l), (107) 

which satisfies 

/ dyi d2/2G^^(yi,2/2)= / dyi dy2(yi + y2)G^^(2/i, yz) = 1- (108) 

Jo Jo Jo Jo 

Consequently, the structure function of the pion J-'j^ix), that is related to the single- valon distribution G^{y) 

T^{x)<x f dyGliy) (109) 



behaves at large a; as (1 — a;) ^, in agreement with the counting rule (1 — x)^'' ^, where r is the number of residual 
quarks (r = 1 for pion and r = 2 for proton). 
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Using Eqs. (ITM) and dlO^) in (WM) yields 



By defining 



1 pl-xi 



JO 



we then have for the moments [see Eq. ( |62|) 



dxi / dx2FTr{xi,X2)S{xi + X2 - x). 



H'^{x)=x^H^{x), 



(110) 



(111) 



pl nl — xi 

H'^{n) = dxi dx2F.^{xi,X2){xi + X2)" 



Jo 



, , n\\n2\ 



(112) 



whereupon Eqs. (104) and (105) can be used for the production of tt"*" and 7r~, respectively. The inverse transform 
can be done as before, using Eqs. ( p8| ) and (|9l|). 

Having formulated the procedure to calculate the pion distribution in the projectile fragmentation region, we must 
confront one final issue on the role of the gluons. Although the gluons carry roughly half the momentum of a proton, 
no glueball has ever been seen. They therefore hadronize by converting to qq pairs, which subsequently form pions. 
We take them into account by enhancing the sea to saturate the momentum sum rule That is, for the purpose 
of pion production we revise the normalization of the quark distribution L{z) such that qq in the sea carry all the 
momentum of the incident proton apart from the momenta of the valence quarks, leaving nothing for the gluons. The 
average momentum fraction carried by the valence quark in a valon is 



^>val 



dzKNs{z) 



where Eq. (^ has been used. If we denote the saturated sea distribution by 

ii(z) =4(1-^)5, 



(113) 



(114) 



where only the constant factor £1 has changed from Eq. (pO|), then each sea quark carries on average a momentum 
fraction of ^1/6. That is to be identified with (1 — (z)yg^]^)/2/, where / is the number of flavors. From Eq. (p^ ) we 
have {z)^g^\ — 0.52. Setting / = 3, we get 



(115) 



With this value in Li[z), which is then used in Eqs. ( ^l| ) - ( |55| ) in place of L{z), we obtain the appropriate V and W 
functions that should be used for the calculation of F^± in Eqs. (104) - ( |105|) . 

Evidently, there are no parameters to adjust for the calculation of the pion distributions. The results are shown 
in Fig. 8 for tt"*" and 7r~ separately at i/ = 3.1 and 6.3. For ease of comparison between 7r+ and 7r~ the same curves 
are replotted in Fig. 9, where the charge and v dependences are grouped differently. No data points are included 
because none correspond to proton fragmentation only and at the v considered, as discussed at the beginning of this 
section. Nevertheless, if one compares our results to the data of E910 shown in (l), there is rough agreement, both 
in normalization and in shape. Generally speaking, the difference between tt"*" and 7r~ diminishes as V is increased, 
though not as rapidly as in the E910 data at E^^^ — 12GeV. We are confident that our predictions will agree well 
with the data when they become available, since the pion distributions have always demonstrated the reliability of 
the recombination model [|12l nsl. 



VII. CONCLUSION 



We have formulated the projectile fragmentation problem in pA collisions in the valon model. Despite the non- 
perturbative aspect of the problem, the formulation results in a well-defined procedure of calculating various con- 
tributions to the nucleon and meson distributions in the proton fragmentation region. The nuclear target effects, 
which present the only unknown in the model, are summarized by two parameters k and q. They are determined by 
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fitting p — p and n — n distributions for two values of />. With those parameters fixed, the predictions for the pion 
distributions can be calculated without any other adjustable parameters. 

The results of our calculations have shown that the NA49 data can be well described by the valon model. The 
normalization of the nucleon spectra turns out to be correct without any freedom for adjustment. The shapes of the 
distributions are also acceptably reproduced. The inferred value of k that gives a quantitative measure of momentum 
degradation can be translated to a degradation length A in the form e"^^-^ for the degraded momentum fraction with 
A ~ 17 fm. 

While a value for A can serve as a succinct numerical summary of the stopping effect of the nucleus, we admit 
that it cannot be inferred directly from the pA collision data without a detailed analysis in the framework of the 
valon model. This aspect of the problem is worthy of further attention in the hope that a degradation length can be 
extracted by an appropriate model-independent analysis of the data on nucleon inclusive cross section. 

It should be noted that since the valon model does not make explicit use of Regge exchanges, it is not capable of 
predicting the energy dependence. It is a model that should be applied only at asymptotic energies where scaling 
behavior prevails. To compensate for that drawback, it makes possible the s-channels approach to the calculation of 
the fragment distributions. 

There are obvious directions into which this work can be extended. One is to incorporate strangeness and study 
the distributions of hyperons and kaons. Another is to generalize from pA to AA collisions. From the properties of 
degradation that this work has revealed, one is better positioned to assess the extent to which nuclear matter can be 
compressed in AA collisions. Furthermore, with some knowledge about strangeness production in pA collisions, one 
can determine for AA collisions how much strangeness enhancement is normal and how much anomalous. 

Although the valon model represents an approach to multiparticle production that appears to be orthogonal to 
most other approaches based on strings, it should be recognized that the s-channel and t-channel approaches are 
complementary, not contradictory. One may be able to identify diagrams in Figs. 3 and 4 that correspond to baryon 
junction or diquark breaking. Just as the notion of duality has benefitted the investigation of high-energy processes 
many years ago that found its idealization in the form of Veneziano amplitude, which has perfect s- and ^-channel 
symmetry, so also here in nuclear processes the exploration of complementary descriptions of common as well as 
unusual phenomena can help to elucidate the underlying dynamics responsible for them. 
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APPENDIX A: APPENDIX 



In Sec. Ill we have considered the problem of describing v coUisions with the target nucleons in terms of the number 
of coUisions that each valon experiences. With the ith valon encountering Vi colhsions, the sum vi satisfies the 
bounds given in Eq. (|l6|). In this Appendix we investigate the phenomenological preference for that sum within that 
range. 

Let us define the integer /i by 

/i = i^l + 1^2 + 1^3 (Al) 

so that [i is bounded by 

V <^i<'iv. (A2) 

Since at least one valon must interact with the nucleus, let that valon be i = 1. Let p be the probability that either 
one of the other two valons also interacts. Furthermore, let B^{ii) denote the probability that out of i/ independent 
collisions the target nucleons encounter, /i valonic collisions occur. For i/ = 1, we have i?i(l) = (1 — p)^,i?i(2) — 
2p{l — p) and -Bi(3) = p^, which count the probabilities that valons i ~ 2 and 3 interact in addition to the i = 1 
valon. Generalizing that to v collisions, we have 

B.ifi) = (^j)p' (1 - P)""'^ J=f^-'^, (A3) 

which is a binomial distribution of having ^ — v valonic collisions by the i ^ 2 and 3 valons out of a maximum of 2v 
possible such collisions. We can define a generating function from By{ij): 

j=0 V / 

The above consideration can now be applied to Eq. ( |3^ ) where G"j,(ni, 712, na) is related to G(ni, 71,2,71.3) under the 
assumption of Eq. ([l^), i. e., i> = vi + V2 + 1^3- That assumption is now liberated by Eq. ( |A1[ ) and (A2). The factor 



d{ni) / 3]'^ in Eq. ( [36| ) should thus be replaced by [Si{z)Y, which allows /x to vary between and 37^. In the 
limit p = 0, Si{z) becomes z and we recover the earlier result, if we identify 



3 

i ^ 

z 



IE^K)- (A5) 



Using ( [A4| ) in Eq. (|36|), we finally have, with the help of Eq. (|2 

Gp(7il, 712,713) = G'(?ll, 712, 713)5^(7/1 + 7^2 + l^3)Pii{v) 

■ 1 \ fl +1^2 +1/3 / I I \ I " 

1\ {Vi+V2+V3r. 



3 / Vl\v2^.l''i\ 



\{d{n,Y. (A6) 



The summation over Vi is included in the sum in Eq. (|2^) but without the restriction of (|17|). 

We have used the Gp given in ([A^) to calculate H'{x) for p — p and n — h, as in Sees. IV and V. The only difference 



from before is that we now have one additional parameter p in (A2) to adjust, which controls the number of valonic 
collisions above ^ = u. Our best fit of the data, as in Fig. 5, yields p = 0.05 with the values of k and q being essentially 
the same as in Eq. (p6|). Since p is so small, any contribution from /i different from v can be neglected. Thus it is the 
data that instruct us on the suitable range of values of vi, namely: vi-\- 1^2 + 1^3 is dominantly at its lower bound v. 
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U Q u V. v; U Q d 
(a) (b) 





FIG. 2: (a) V/^, favored quark in a valon, (b) V^^, unfavored quark in a valon, (c) W/^^, two favored quarks in a valon, (d) 



Wf^f,, a favored and an unfavored quark in a valon. 




FIG. 4: Eight types of contributions to the quark state udd in a proton 
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FIG. 5: Inclusive distributions for the production of p — p (solid) and n ~ n (dashed) in p-Ph collisions. The data are from 
NA49, reported in jlj. 
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FIG. 6: (a) A dual diagram that represents the process A"*" ^ n + tt^. (b) The square box symbolizes a flavor change from U 
valon to D valon before the quark momenta are determined in the valon model. 




FIG. 7: Five types of contributions to the quark state ud in a proton 
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